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Abstract 

The  problem  of  diffraction  of  a  plane  wave  by  a  circular 
aperture  in  a  plane  screen  has  been  treated  by  Levine  and  Schwinger, 
who  reduced  it  to  a  system  of  infinitely  many  linear  equations. 
This  system  has  been  modified  by  Bouwkamp,  and  Magnus  has  proved 
that  Bouwkamp' s  system  is  equivalent  to  a  Fredholm  integral  equa- 
tion of  the  second  kind.  If  a  =  ka,  where  k  is  the  wave  number 
of  the  incident  wave  and  a  denotes  the  radius  of  the  circular 
aperture,  then  the  coefficients  of  Bowkamp's  sytem  depend  on  a. 
Their  asymptotic  behavior  for  a  ->oo  is  investigated  and  explicit 
expressions  for  the  first  few  terms  are  stated  in  this  report. 
The  integral  equation  equivalent  to  Bouwkamp' s  system  degenerates 
into  an  integral  equation  of  the  first  kind  as  a  ->oo.  The  solu- 
tion of  the  degenerate  integral  equation  can  be  given  explicitly, 
and  a  perturbation  method  yields  some  formulas  for  the  solution 
of  the  original  integral  equation  if  a  is  large. 
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1*  Introduction 

We  consider  here  the  problem  of  the  diffraction  of  a  plane  scalar 
wave  normally  incident  on  a  plane  screen  with  a  circular  aperture.  Let  u  be 
the  solution  of  Zlu  +  Ku  ■  0,  which  represents  the  diffracted  wave,  and  let 
p,  ©  be   polar  coordinates  on  the  screen  such  that  p  ■  0  represents  the  center 
of  the  aperture.  We  assume  that  on  the  screen  u  =  0.  Let  the  radius  of  the 
aperture  be  a  and  let  $  (p)  be  the  value  of  u  in  the  aperture,  0  <  p  <  a. 
BouwkampL-',  using  the  Levine  and  Schwinger'-  ■*    variational  method,  derived  the 
following  result. 

If  $  (p)  is  expanded  in  a  series  of  Legendre  polynomials  of  the  type 


(1)         $  (p)  -  £  bn  P2n+1  U   -  p2/a2^) 


,  n  2n+l 
n=0 


then  the  b  satisfy  the  system  of  infinitely  many  linear  equations 
(2) 


^0  m,n  n    xa    m,o 


where  a  =  Ka,  6   =1  and  6   =  0  if  m  4  0,  and  where 

'     o,o  m,o  '      ' 

•6\2     Rn  +  I)        Rm  +  §)       * 
(3)  d    .  -     (?)      -L      -1-     g!  Ja) 


m,n        W  ,  .  ^ 

*  s  '  n.  m. 


oo 

r 


5a,n* 


V 


^^     J2m+3/2(av)   J2m+3A(av)  dV* 


In  (I4) ,   J  denotes  the  Bessel  function  of  the  first  kind,  and 

(Si  /v2- 1   -   -  iyi-v2',  >i-v2'  >  0 

for  0  <  v  <  1.  The  diffracted  amplitude  in  the  forward  direction  is  given  by 


(6) 


-  2  - 


Instead  of  using  Bouwkamp's  coefficients  b  for  the  expansion  of  the 


field  in  the  aperture,  we  introduce  the  quantities 

n  rw§) 

(7)  s  =  (-l)n  S_   b   . 

v  n  •  _i      n 

n. 


oo 


/ § — 21 

Let  y  »  yi  -  p  /a  ,  so  that 

(8) 

Then  Magnus  has  shown ^ -*  that  ^3  (y)  satisfies  the  following  integral  equation: 


n=(J 


2a 


(9)  2HL     =     fe(y)    +^t 


-iay 

y* 


; 


G(x,yja)   ^(x)   dx, 


0  <  y  <  1 


The  kernel  G(x,y)   is   given  by 


(10)     G(x,y,a) 


x+y 


r 


dr  , 


) 


|x-y| 


where  Jn  is  a  Bessel  function  and  H.  is  a  Strove  function  and  a  is  a  real  parameter. 

From  a  knowledge  of  &(y)   we  can  (using  (8),  (7)  and  (1))  obtain  |  (p), 
the  value  of  u  in  the  aperture. 


2.  The  field  in  the  aperture  for  a  ^  oo 

Magnus  has  shown  **•*   that  a  solution  of  (9)  exists  for  all  real  positive 
a*  We  wish  to  study  the  asymptotic  solution  of  this  equation  for  large  postive  a. 
Our  first  task  is  to  determine  how  the  solution  <§(xja)  behaves  for  a  -*oo.  To 
do  this  we  will  need  to  know 
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lim      G(x,yja) 
a  ->oo 


5    G(x,y;oo) 


for  0  <  x,  y  <  1  and  for  x  /  y.     From  (lO) 
G(x,yj  oo) 


lira 
a  — »  oo 


n 


x+y 

f  Jx(ar)    +  iR^ar) 

_ 

J  N-y| 


lim 

i  — ►  oo 


x+y 


|x-y| 


Jx(ar) 
r 


dT     »      lim      -  J  JJor) 
a  — *•  oo 


dr   . 


x+y 


; 


|x-y| 


dr 
r 


-flog 


2^2      lim      JUar)   , 

x-y    _     .  „    \L  ' 


a  ->oo 


where    |x-y|  <     Z(a)  <  |x+y|.     Since  x  /  y,   J,(ar)  -^Oasa  -*oo.     Likewise 


lim 
a  -»-oo 


x+y 


y 


|x-y| 


r 


dr 


-  rlQg 


x+y 


x-y 


lim      ^(ar*)  , 


a  ->oo 


where   |x-y|  ^    T  (a)  <  |x+y|.     Since  x  /  y, 


lim       H^ar*)     -  j1       1       "    f 

a  -»oo  I   (J)    I   V£) 


Hence  if  x  /  y 

(10a)  G(x,yj  oo)     -    -  §  log    |£ 


for       0  £  x,     y    Si 


Let  S>  (x)  denote  the  limit  of  <S(x$a)  as  a  -►oo.  We  can  now  show 
Theorem  1»  For  a  -»  oo,  the  solution  of  (9)  tends  towards 

g(x)  .  A.     __2 —  . 

^    y5T? 


-1»- 


Proof;  Divide  (9)  by  a  and  let  a->oo.  We  assume  that  lim   JgiZi£J  =  o,  and  that 


a-*oo 


lim      G(x,y$a)  S(x,a)dx  ■    G(x,yjoo)  S  (x)  dx  | 
a  -*-oo  I  i 

the  validity  of  these  assumptions  will  be  justified  later.  Then, taking  into 

account  (10a),  we  have  the  following  integral  equation  of  the  first  kind  for  the 

determination  of  <§  (x)  : 

o 


(11) 


A 


log 


sz  £  (x)  dxj 

x-y|  o     ' 


0  <  y  <  1 


The  easiest  way  to  obtain  the  solution  of  (ll)  is  to  proceed  from  the  known  solu- 

DG. 


tion  of  the  following  integral  equation 

+1 

1 


(12) 


'hfi? 


log|x-y|  dx  ; 


-1  <  7  <  1 


This  may  be  transformed  as  follows: 

? 

x 
7 


1 
n 


J-i   >C7 


log|x-y|dx  -  - 


5  '.     £7 

1 
1      / 


log|x+y|dx  -  - 


o    -pL-x 

1 

r      x 


>G? 


log 


x+y 


x-y 


n 


dx     • 


>C? 


log|x-y|dx 


log | x-y |  dx 


Comparison  of  this  result  with  equation  (ll)  yields  the  desired  result 


Six)    .  -i-  _5_ 


-5- 


For  this  limiting  case  we  have  the  first  transmission  coefficient  A_ : 


*1 


-^b 
3  o 


s 

ia    o 


-  ia 


1 

r 


T 


nf )     3  nf ) 


•  3 


P1(x)  £0(x)dx 


-2ia 


7o 


1   X   . 

x  —  ■     dx 


ia 

T  * 


We  now  justify  our  assumption  that  letting  a  -»oo  while  keeping  the 
radius  of  the  aperture  constant  is  equivalent  to  keeping  a  constant  while  let- 
ting the  radius  of  the  aperture  approach  oo.  If  a  is  constant  and  the  radius 
of  the  aperture  tends  to  infinity  ve  know  that  the  wave  goes  through  undisturbed. 
Hence  if  &>   (x)  makes  ^  (p)  constant  in  the  aperture  we  feel  justified  in  our 
assumption.  We  will  now  show  that  this  is  indeed  the  case. 

We  have  from  (8) 

oo 


(13) 


g0(x) 


yn  £7>  n=0 


£  **   P2n+l(x)   ' 


Using  the  orthogonality  property  of  the  Legendre  polynomials  we  get 


n 


o  ^7 


r 


kn  +  3 
Un  +  3 

u  yff 

/  o 

hZ-Zl.  2nl 
U  /5T  Vn, 


P2n+l(x)  <* 


P2  +i^cos  e^   cos  e  de 


'2n\     /2n+2N 


n+1 


,-Un-2 


See   [U],  p.  $1, 
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kn+3     y~  R2n+l)    Qg2g) 

2  n!ni(n+l)iCn+l):  2Un+* 

Un+3  Rn*|)n:  22n+1   Rn  +  f)(n+l)j  22n+3        t 

2  "       (n.')2(n+l)i(n+l):   2hn+2  2  /n     .   2  /JT 

Un.3         Rn  +  |)    Rn+§) 


2  /n  n!   (n+1)! 

From  (7)  we  get 


Rn+i)    Rn  +  |) 


n         Rn  +  |)       2  ytf  nJ   (n+1): 

.  (-D-  k-1  r(D :  ij  . 

2  y?T       Rn  +  2) 


Let  us  now  determine  the  C  in  the  relation 

n 


1    'A    °n  W*> 

n=0 

• 

get 

C     ■ 
n 

1 
WH     f1'   P2n+l(x)   to     = 

fell  (-1)" 

1 
2" 

3          2n+l 
2   ***  " 2    ' 

(n+1)  .» 

X 

/  ^  n    Un  +  3 

Rn  +  |) 

2  y5     Rn  +  2) 

Hence  C  =  b  and  J  (p)  is  constant  in  the  aperture. 
n   n     * 


Here  we  have  used  the  formula  R2z)  =  -i—  22z  Rz)  Rz  +  i)  . 

2  /n 
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There  exists  a  general  formula  which  permits  us  to  express  the  field 
($  (p)  in  the  aperture  in  terms  of  the  solution  S(y)  of  the  integral  equation 
(!?)•  The  connection  between  the  two  functions  is  given  by  the  following: 

Lemma  1:  Let  G(y)  be  given  by  (8)  and  let 


(x)  =  £2  (-Dn 


t£o       rw£) 


7"      n  r2n+l 
2 


P^^U)    . 


Then 


|  (x) 


1 

r 


n  U,y)    £(y)  dy  , 


J 


where 

n(x,y) 


1 

i  A  J 


; 


and 


W(x,y,t)     -    ^ 


n 

r 


r(x,y,it)    -  W(x,y,-it) 


(1  -  t  )     do 


dt 


>C? 


{^[xy.d-x^^Cl.y^^cos  .]   +  t2} 


372**   * 


H 


Proof;  According  to  Watson LJ, 

oo 


w(x,y,t)   =  £]  (n  + |>  pnw  pn(y) tn  • 


n=0 


Therefore, 


1      Jw(x,y^  -  W(x,y,-it)  I    =    £]     <2n  +  ¥   P2n+1^>   P2n+1^>    W  £"*  • 


2T 


n-0 


Since 


.2n+l       dt 


y& 


we  see  that 


n(,)7)  -  s  (-Dn  -Caa  P    w  Wy>. 

I  ln  +  »; 


This  and  the  orthogonality  relations  for  the  Legendre  functions  prove  Lemma  1, 
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3»  Asymptotic  expansion  of  the  diffraction  integrals 

In  this  section  we  study  the  asymptotic  expansion  of  the  diffraction 
integrals  (h)*     (U)  may  be  written 

g^n(a)=-i         J2^2_  J^/gUvJJ^/^aOdv  ♦     I      ^1  J2ra+3/2(av)J2n+3/2(av)dv. 


J 


v 
o 


1         V 


We  treat  the  ±'eal  and  imaginary  parts  separately.     Consider 

Ira*     C)  ? 

-Imeni,n(a)     ' 


JT/Z 

'   g         J  (av)    J  (av)  dv 

v 
o 


3  3  + 

where  n.»2ra  +  »,  v»2n  +  *.     We  have ' 

(-l)n(°^+v+2n   n^+v+l+2n) 
J  (ov)J  (ov)     =  Z2 


n-0      nJRu+l+n)    Rv+l+n)    Ru+v+l+n) 


n=0      ni       V2/  2  yn   Rn+l+n)    Rv+l+n)    Ru+v+l+n) 


Substituting  this  sura  in  our  integral  and  interchanging  summation  and  integration 
we  obtain 

^>n  yE    ^0      n!         Rn+l+n)    Rv+l+n)    P  (n+v+l+n)  I 

Evaluating  the  integral  on  the  right,  we  have 

yC7  ^+v"2+2ndv  -    I    (sine)  V2^    yC0s29de  »       gp7gg ;  f    ■ 


tsee    (6j,  p.  Iii7» 
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So,  finall7, 

-L,  g*     (a)  -  V-  7^    i^J LJCH LJ_iZL_Z  «  °  V]    iz2-L  f(n), 

^n  h       1^0       nl  P(^l+n)Rv+l+n)H^+v+l+n)  *»      Vo      nl 

where 

P  (kt+l+n)  r  (v+l+n)  P  (u+v+l+n) 
To  obtain  the  asymptotic  expansion  of  a  series  of  the  form 

fl     ^4^    f(n)     , 

n-0        n; 

for  large  positive  o  we  use  a  method  due  to  E.M.  Wright1-  J. 

Using  Wright's  notation  we  now  state  that  the  following  asymptotic 
expansion  holds: 

-  Im  g*  n(a)  ~  1(2^)  +  I(Z2)  +  J(a2)  , 

where  the  meaning  of  these  symbols  is  defined  for  our  case  by  (for  details  see  [7]])  : 


I(x)  ■  x  e  < 


f  M-l 


m=0 


*.  -  4  £  >„  f**»*>  •  0(O  . 

In  the  above  equations  M  denotes  any  positive  integer;  N  is  an  integer  which  may 
be  taken  as  large  as  we  pleasej  and  6  is  an  arbitrary  positive  number,  and  finally 
L  is  an  arbitrary  integer.  If 

P  (p.*  a. t)  T  (P2+  agt) 

f(t) 


P  (^  f^t)  P  (,y  p2t)  P  (iy  P3t) 
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has  a  pole  of  order  s  at  the  point  -(^+6  )/a  ,  then  s  P    y    r'/**r  ±s   the 


r  '  r' 


residue  of  I_l(-t)f(t)y  at  this  point;  if  a  >  1,  then  P  .  is  a  polynomial  in 

r,x 

log  y  of  degree  s  -  1.  We  know  that 


r(pto)  -  e-°t(at)P+at-1^  < 


M-l 


-Ma 


S   *   t-  ♦  0(0 


m-0 


ra 


>   • 


Using  this  fact  we  arrive  at  the  following  inequality: 


f(t) 


M-l 


u  n t+i)     »*  n  2t-e+m+i) 


K 


P(2t-€+M+l) 


where  K  is  independent  of  t. 

This  inequality  serves  to  determine  the  number  A  ,  A,, 
In  particular,  if  we  let  M  ■  1  and  let  t  -»  oo,  we  get 

ao=  22M*2n+6/y2«    . 


. .  uniquely. 


Using  this  procedure  we  can  define  further  A. 's  recursively, 
linally,  we  need  the  following: 


2,   q  -  3,  o^ 


P!  -  P2  -  P3  -  1 


Pi "  — r^  » 


b   =  ^T1 

P2   *T^ 


1^  -   n+l,    n2  «  V+l,     U^  -  ti+V+1 


*  -4>r-£.v¥ 


r-1 
^  -  2cd, 


-2m-2n  -  ■*■ 


Z2  =  -2ai  • 


Now  f(t)  has  a  simple  pole  at  t  -  -  8-j/ol  -  -(n+v-l)/2  .  All  its  other  poles  are 
of  order  2  and  occur  at 
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l£j   .  .  (ustpm       tBt      t .  1A  ...  . 


V     i  "5 

The  residue  of  PC-t^t)^  at 


t r- 


is 


r^u+v-lj  -(m-+v-1)/2 
r(ii!vi3)  ^^j  r(v^3}   • 

The  residue  of  Q(t)  =  PC-tJf^y*  at 

is  obtained  as  follows.  We  first  obtain  the  Laurent  expansion  of 

H?!*  c^t)  =  Rp^  t) ,   at  t  -  -(p.^  J? )  o 

Let  Z"  p..+  t.  Then  Z.-  -  i   when  t  ■  -(£-,+  i),  and  we  now  need  the  Laurent  ex- 
pansion of  Rz_)  about  Z--  -#.  We  obtain  (cf.  [8],  p.  16) 

Next  we  need  the  Laurent  expansion  of 

P(p2+  ogt)  -  Rp^  1  +  t)   at  t  -  -(p-^JO  . 

Let  Z2»  p  +  1  +  t.  Then  Z2=  1  -  £  when  t  »  -(p,+  J?),  and  now  we  seek 
the  Laurent  expansion  of  Rz  )  about  Zg-  1  -  J?  -  -( J?  -  l)j  we  obtain 


Rzj  -  (-1}    Hz9*S.  -1)' 
2     (J?  -Di 


*  +  T(J!)  + ...  | 
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We  now  obtain  the  expansion  of 


R-t)/ 


r<iy  t)niy  t)r(iy  t) 


p(t) 


at  the  regular  point  t  ■  -(p. +  $  ): 

P(t)     =    P(l(p1+J?))     +P,^(p1+I))    (t  +  (p1+i?j) 

[V I  )y"(pl+  * ) 


Now 


and 


_d_ 

dt 


R-t) 


Ry^H^  t)Riy  t) 


-(P1+  f ) 


t  -  -(p^n 


R-t) 


R^*  t)Rix2+  t)Rty  t) 


-0-,+  « ) 
log  y  y      J- 


t  =  -(03+i) 


(t  +  px+  J?  ) 


R^)     -    iiL    j(t  +  P1+  ft)"1  +  ¥(l  +  1)  +  ...  ]  , 

Rzj    -    W  "1t    j(t  ♦  p,+  *)"1  +H!L)  *  ...  X, 
2  (5-D!    \         1  J 


P(t)     -    y 


-(p1+ *  ) 


Rpx*  < ) 


Rv  pr n n^-  Pj,- 8) Ry  tfl > 


R-t) 


dt/     3 

TT  Rm  t) 

i-l 


,  3  r  c-t) 

wp^  &rv*> 


log  y 
t=-(p1+^) 


(t+p1+«)  ^ 


-(Pi*  * ) 


A  +  (B  +  C  log  y)   (t  +  px+  4 ) 
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where  the  meaning  of  A,  B,  C  is  obvious.  Finally 


^      2       UK  a  -d: 


i      +  y(j?)  + *o?+  i)  ,  _ 


Hence  the  residue  of  Q(t)   -    R-^fCtJy*    at  t  -  -(px+  Q  )  is 

^"1)y A  (Y(tf)   +  T(*  +  l)j       +  (B  +  Clogy)         . 


With  these  calculations  before  us  we  write  the  asymptotic  expansion  of 


-  Im 


\%n 


(a)        -   T\ 


00        (^2)n 


n-0        ni 


f(n) 


for  large  positive  a.  It  is  given  by 


1(2^)  +  I(Z2)  +  J(a2)   -  ££L    J(2ai)ee2ailr;Am(2ai)-m  +  Q  (W)-*) 


♦  (-2ai)V2ai 


p.+v 


Tm 

Z:  Am(-2ai)- 
l_m»0 


♦  0  (a"M) 


""ir  i  njj^  r^9  rt^)  3^ 

1*1^37?^  [+(B  +  2cloga)    J'  • 


Y(Z)  is  the  logarithmic  derivative  of  the  gamma  function: 


rm  -  4  log  r  w  ■  £gj.   . 


-  Ik  - 

We  now  consider  the  real  part 


*2p  J^Cav)  Jv(av)dv  -  Re  Jg^a) 


J 


We  remark  first  that  if  only  a  few  terms  of  the  asymptotic  expansion  of  the  real 
part  are  desired,  we  can  obtain  these  by  substituting  for  J  (civ)  and  Jv(av)  the 
first  few  terms  of  their  known  asymptotic  expansions  (in  fact  these  expansions 
consist  of  a  finite  number  of  terms).  In  what  follows  we  give  a  slight  modifi- 
cation of  this  procedure  in  which  we  will  again  make  use  of  Wright's  methods. 
VJe  have 

J  (ay) J  (av)     =     (av)  72      £2JlL      -LA-? Z-Lkj -Z_      . 

V-         v  ?Z         rA         n!  P  (n+l+n)  P  (v+l+n)  T  (n+v+l+n) 

Since  the  range  of  integration  is  from  1  to  infinity  we  know  that  av  is  large  if 
a  is  large.  We  obtain  now  the  asymptotic  expansion  of  this  product  for  large  av. 
Let 

r  (ji+i+t)  r  (v+i+t)  r  (n+v+m) 

The  following  will  be  needed: 

p     «     2,       q     -     3,       Oj^  =     a2  »     c^  =     p2  -     p3  -     1  J 

Pi  '    ■nr  •  P2  "    ^    J 

u^  =     M.+1,       jx2  =  v+1,       |Xj  -     H+V+l  J 

©     -     II  V  i    ^r  +    I  (q-p)     "     "  ^  "  2n  "  U  5 
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g(t) 


M-l 

22 


A 


m 


It  r(t+1)     m°°   P(2t-Gtm+1) 


P(2t-©+M+l) 


where  K  is  independent  of  t.  This  last  inequality  makes  possible  the  calculation 
of  the  A±.  In  particular  A  -  1/  yn  22m+2n+il.  Also 

Z.  -  2avi  ,     Zg  -  -  2avi  . 

g(t)  has  simple  poles  only  and  these  occur  at  t  ■  -(p2+^)>  $   =  0>3-»2,..., 
(m+n+2). 

The  residue  of  R-tOgMy*  at  t  -  -  (pg+  J?)  is 

j  .(E2|a+j) 


Thus  the  asymptotic  expansion  of  J  (av)J  (av)  for  av  large  is 


(ov) 


u+v 


yn 


J(2avi)Q  e2avi  22   A  (2avi)-m  ♦  0({a*)^) 

Dl=0  ^-     ^ 


+  (-2avi)6  e"2aviS  Am(-2avi)Ha  ♦  O^v)"^ 


_.  M 

m=0 

m+n+2  (-1yi         RJ2g&*nn!-n 


rU-V+1 


Vi+v+1 


p3    «     p(i^.fi)n^p-nn:i2p-n(aT) 


.u+v+l+2j2 


We  now  multiply  by   '  T  .   and  integrate  from  1  to  infinity.  Then  we  get  es- 

v 
sentially  two  different  types  of  integrals  to  evaulate.  The  first  type  is  of  the 

form 


J 


1 
f 


v-1 


T2T 


dv 


.28 


>C7 


dt 


?       2o   2        rh  +  |>ri  +  |> 

(cos  e)*xsin  ©d©    - 


2T(  I   +2) 
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The  second  type  is  more  difficult  to  evaluate.  It  is  of  the  form 


oo 

r 


£H 


We  let  v-1  «*  t  and  we  get 

_2cd 


(av)P 


2avi   . 
e  dv  , 


p  an  integer. 


oo 

r 


2SpZ    e-(-2ai)t    dt     . 


A 


(t+1) 


Sp 


An  application  of  Watson's  lemma  (see  J6],  p«  236)  will  give  the  asymptotic  ex- 
pansion of  this  last  integral.  We  see  then  that,  after  some  tedious  calculations, 
we  can  obtain  the  complete  asymptotic  expansion  of  the  diffraction  integrals 
e   (a)»  For  the  first terms  of  these  expansions,  we  have: 

Theorem  2.  Let  e   (a)  be  the  functions  defined  by  equation  (U), 


,n 


to  terms  of  the  order  a 
vm+n 


and  let  a  -»oo.       Then  the  asymptotic  expansion  of  the  g_     (a)  up 
e  ore 

(-1)1 


Re 


Im 


**m,n 


(a) 


TioT 


is  given  by 


cos  (2a  +  4~) 


jr 


^n 


(a) 


P(m+n+2) 


k  r(™+n+3)  r(m-n+3/2)Rn-m+3/2) 


+  i^ 


n+m 


2n 


n+m 


a     +     (l-2y)(-l)*™  1 
lin"  a 


B(t)     1     ,      (-1) 

u       a  o  ^ 


m+n+l      cos  (2a  -  4~) 


2  yn 


where  y  denotes  Euler's  constant  and  inhere 


W 


-=*—    +  ...  , 


B(t) 


dl  i 


R-t) 


n2m+V2+t)  f[2n+S/2+t)  P^+^+h+t) 
evaluated  at  t  «  -  (m+n+2). 
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U»     An  approximation  formula  for  large  a 

In  this  section  we  study  the  next  term  in  the  asymptotic  expansion  for 
a  -»co  of  the  solution  G(x)  of  our  integral  equation 


(9) 


-*S  -  g(y)  +~  f  G(x,yja)g(x)  dx  . 


To  obtain  this  information,  we  integrate  (9) •     We  get 

dy     ■ 


z 

.is: 


z  z       1 

f      _  2a      C        r  6 

S(y)dy  +  ~  G(x,yja)  &(x)dxdy 


J 


J 


o    J 


or 


(15a) 


where 


and 


iaz 

2  ya 


E(z)  +  ^  K(z,x)  e'(x)  dx  , 


z 


K(z,x) 


0  <  z  <  1  ,  E(z) 


G(y>x)dy     =     -  £ 


S(y)  dy 


f  jS? 


iaw 


[Z  iawx  iaw(x+z)    _  1 


+  sgn(z-x)(e 


iaw  z-x 


-  1)1  dw     . 


For  later  use  we  note 


dK 


i 
2 


1 
f 


TJ-w       L  iawx  m       iaw(x+z)_       iaw|x-z| 


} 


dw  • 


An  integration  by  parts  of  (l£a)  gives 


(15b) 
where 


xaz 


E(z)  + 


2a 


K(z,l)E(l) 


z 


2  /n 

K(z,<r)  -    G(a,y;a)  dy  -  0 


-  J0iE(x)<*j  ' 
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since 


G(<7,y;a) 


J0^'    j^or)  +  i^Car) 


dr    2   o   . 


J\<y-7\ 


Equation  (l£b)  is  exact.  The  limiting  solution  as  a  -»oo  will  therefore  be 
given  by 


Eo(z)  - 


z 

a.  g  (x)dx  =  -  A.  >C7  +  A 


We  wish  now  to  obtain  some  information  about  the  next  term  E-(z)  in 

the  asymptotic  expansion  of  E(z).  Towards  this  end  we  now  make  the  following 

approximations  in  equation  (lf?b).  We  write  E(z)  =  E  (z)  +  E^(z)  +  ...  and 

substitute  this  into  (l5b)  dropping  terras  of  higher  order.  We  therefore  replace 

E(z)  by  E  (z),  E(l)  by  E  (l)  and  E(x)  (which  occurs  in  the  integrand)  by 
o  o 

E  (x)  +  E  (x).  We  get 

-  1 
o 


iaz 


2a 


(16) 


2  yn 

2 


Eo(z)  -  £  K(z,l)  ±-    ♦ 


|£  E  (x)dx 
3x   ox 


iaz 

2 


z     t-,  /  \    2a  „/  n\    2a     3K 


dx 


2a 

o 


i£? 


dx 


2a 
n 


1 

r 


1  y*>  *■ 


y 


This  is  an  integral  equation  of  the  first  kind  for  the  determination  of  E^(x) , 
since  the  left-hand  side  contains  known  quantities. 

We  now  show  that  after  division  by  a  the  left-hand  side  of  (16)  vanishes 
to  a  higher  order  than  a"5'  as  a  ->oo.  For  this  purpose  we  need  the  following 
calculations: 


-  1°  - 


2a 


1)     mm£L.  K(z,l) 


2eiaw  m  eiaw(l+z)  _  eiaw(l-z)  . 


'} 


1 

2a* 


i 
-2a 


r 


£7 


m 

Jo 

1      , 


'2  cos  aw  -  cos  aw(l+z)  -  cos  aw(l-z) I     dw 


2  sin  aw  -  sin  aw(l+z)  -  sin  aw(l-z) I    dw 


1 


i 

2a- 


; 


w 
o 


/"./■>     •    /       .  ccwz\    .     awz\      /"0     •    /■         awzx    .     awzNl 
(+2  sin(aw  +  -s-)sin  -r-l  -(2  sin(aw-  -^-)sin  -75-H 


1 

r 


dw 


-w 


/ 


/%         /       .  awzv 
/-2  cos(aw  +  — x— )sm 


-r-;sm  -*— 


-)  +  (2  cos(aw-—s— )sm 


awz\_.  .  awz> 

"TV 


awz  .      awz 

dw       '    p —      sin  -p--    cos  aw  sin 

w 


22    K   J^ 
"IS. 

JO   1 

I  w 

y° 

-  Ui     f       .   2     awz     .  >Cw         . 

+  w-  sin      -7j~  sin  aw    -*-=tj —      dw 

Jo  W 


"T" 


>l-w         .     awz     .     _      .    /awzx 
dw      '   n        sin  — x-  sin  aw  sin(— x— ) 


^ 


1 


2     awz 


sin      — jt-  cos  aw      '    g —      dw     • 


2£E 


Therefore 


2a 


K(z,l) 


n  yn 


.   2  awz  aL-w    j         j 

sin    — w-  cos  aw     '    a       dw  +  i 

w 


£■ 


.   2  awz 

in    -tj—  sin  aw 


£7 


dw 
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2) 


we  need  the  value  of  — rrr    T5T  dx'     Calculations  2)  and  3)  give  its  value: 

"n~  (~  ?M+  — )    **    %F  '      C2  cos  awc  "  cos  ^z+x)  ~  cos  aw(x-z))  d: 
Vyn  '  Jo     W    Jo  ^  y 


la 

7* 


w    ]  V    a 


yo 

1 


sin  awx  _  sin  aw(z+x)   sin  aw(x-z) 

aw  aw  +aw 


.1 
o 


7o 


dw   ^"g  ■  J   2  sin  aw  -  sin  aw(z+l)  -  sin  aw(l-z) 
w 


1 
r 


7* 


dw 


-w 


2  sin  aw  -  2  sin  aw  cos  awz 


w 


21 

^7? 


; 


dw 


-w 


w 


1  -  cos  awz^    sin  aw 


Id 

7F 


yC7 


.2     awz      . 
sin  aw  sin      —tt-    dw  $ 


3) 


t  (-  §>(♦#>  J* -^  £  (! 


sin  awx  -  sin  aw(z+x)  -  sin  aw|z-x|  )  dx 


-D 


. . .   dw  .J..r~ I, ,  J- 


n  yn 


w    ^       aw 


cos  aw 


(z+x) 


aw 


cos  aw(z-x) 


-  aw 


z    cos  aw(x-z) 

aw 
o 
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1 


1  —      (  1 

r  -JL     2       -2  cos  aw  +  2  +  cos  aw(z+l)   -  cos  awz 

dw  Co.  J  I 


*      f»-rf? 


n  y?r 


y 


w 


-1  +  cos  awz  +  cos  aw(l-z)  -  1 


-  2  cos  aw  +  2  cos  aw  cos  awz 


J-  fL-rf? 

n  7n  w 


.    1  -  cos  awzL    cos  aw 


n  yz 


dw 


jg 


.    2     awz 
sxn       -*—    cos  aw     • 


We  note  that  l)   +  2)   +  3)   s  0.     Consequently  equation  (16)   simplies  to 


(16a) 


Vjj  iz  n  (    v 

"T        +    2a    Eo(z)     + 


1 

y5r 


dK   7^? 


; 


3x 


dx 


y 


I    Vx>*     • 


We  consider  now  the  integral  in  the  left-hand  side  of  equation  (16a) : 


1 


} 


§>c? 


dx 


2  ya 


w 


} 


dx 


y 


>^{o 


2cos  awx  -  cos  aw(z+x)   - 


-  cos  aw(z-x))  +  i.\2   sin  awx  -  sin  aw(z+x)  -  sin  awlz-xQ 


>  . 


We  wish  to  know  the  value  of  the  imaginary  part  of  this  last  expression  as 
a  ->  oo.  We  have 
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i 
2  yn 


1  /— ^      1 

dw    ■  * ""*    ■  )a-x        (2  cos  avoc  -  cos  aw(z+x)  -  cos  aw(z-x)j 


dx 


J 


o  J° 

2  /n" 


dv      ^  "■■■     U  sin 
w 


2     awz 


/■ 


cos  awx 


>C? 


dx 


7 


-  iyn 


-  i  Vn 


>*  j*2 


y 


a 

r 


sin 


aw 


2  tz 


7T    "1 


T  /     v      .2     awz 
J,  (aw;  sin      -*— 


Jn(t) 


*F  >^F-  • 


; 


The  limit  as  a  -*  oo  of  this  last  expression  is 


-  i  Vn 


and 


J 


P       2  tz  ¥t} 

sin  -K-    -=75 —  dt  •=  R(z)  , 

t 


R 


'w  -  -^ 


sin  tz  — r —  dt 


-4^5 


(cf.  [iQ,  p.  36).  Hence 


R  (z) ^ — 


and  this  cancels  the  first  terra  of  the  left-hand  side  of  equation  (16a). 

We  now  differentiate  equation  (16a)  with  respect  to  a,  and  evaluate 

We  have,  considering  each  term  separately, 


we  now  en: 


-  23  - 


A(a)     ■     2i  J 


7l^w      A-x    dwdx  +  i   j      I  x  sin  awx  ^1-w    7l-x    dwdx  I  $ 

/ojo 


¥a)   -njliSiffa,        0 


=£     f      J>>     £7 


dx 


-ni 
a 


-ni 


J    (ODC) 

o 


£?N 


+    i     [     J  (ax)     — *—      dx 


.i.i  r 

a      a 

Jo 


J  (ax) 


£? 


dx   \.      I 


A2(o) 


-  2  — 

*  2a 


"   H^ax)  >C?  dx     -    =2.    J      J     |  }C?  dx  -     j     H^(ax)  >C7  dx 
o  I    Jo  Jo 


£  n      n 
a  IT      a 


H  (ax) 


>C? 


a 


; 


H  (ax)     ~ —     dx 


£7 


B(e)     .-iff    (x«)eiau(:"2)     >C7    )C7    to* 


-  i 


(x+z) 


0  Jo 

1  1 


(cos  aw(x+z)   +  i  sin  aw(x+z)  L    )l-w      /l-x       dxdw  $ 


}y& 


h 


-  i 


(x+z)  cos  aw(x+z)     )4-w      }a-x      dxdw 
f     J  (a(x+z)) 


o  y  o 


"i     f      "1 

2  L  — - 


■  + 


ni 
2a 


Vo(a(x+z)) 


)C? 


)C? 


t    I  Jo(a(x+z))  —2—  dxl  ^-  2ij  J0(az)j 


£? 


-  2h  - 


1     1 

r  r 


B2(a) 


>7 


n         n 


,7 — 5"     -/""""J  n      I     H.  (a(x+z)) 

(x+z)   sin  aw(x+z)     7j-~w       A~x     dwdx     =     +  7 


£? 


dx 


l-fc? 


dx    - 


H^(a(x+z)     /C?    dx 


'h  (a(x+z)) 


>G7 


1 


1 

+     — 
a 
0  J 


H  (a(x+z))    — ^—    dx 

0         tC? 


n  n 

+  i —  +  — \ 

U      2a 
1    1 


H  (az)    • 


C(a) 


|x-z|   l 

'  J° 

1 


cos  aw  x-z 


-  l 


r     ^(alx-zj) 


OjU)    =   - 


ni 


25 


"i 

2s 


+  i  sin  aw|x-z|l    Vl-w      )4-x      dwdx 
yC?  dx  +  1     j       -L__ yC?    dx  j 

^(aCz-x))     /1-x     dx     +  ^(aCx-z))     yi-x       dx 

/(a(z-x))     yCx    dx     +  /(a(x-z))   ?l-x         dx 


'j  (a(z-x)) 
o^ 


-  a 


>G? 


a     I       o 
Jo 


J  (a(z-x))       -  dx 


Jo(a(x-z)) 


ni 
2a" 


.  j£7.    ♦  ii±  .  j£7 


>G? 


?C7 


?•  i 


dx 


^ 
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C2(o) 


55 


n         n 

Tia  "  2a 


\    ?C?dx-     fV    (a|x-z|)     )£?    dx 
o  Jo 

}  , 


(a|x-z|)     7C?     dx 


n         n 
Ua  "  2a 


JL  _    n 
Ua  "  2"a 


z  J-  ^_^ 

H'(a(z-x))     )Cx    dx     +  H^(a(x-s))     )0-x       dx 

o 

H  fa(SHx)) 


-  a 


£?1    -  i    h0c.(«)) 


£7 


dx 


\(«(m»     ^ 


J 


H  (a(x-z))     -§—    dx 
£7 


1  Ho(az) 

Ha  "  5a    1        a 


Combining,  we  find  that 


JL 
2/ff 


_  (A  +  B  +  c)   -  4  -  4  >C?  ♦  R 

n     L  2a2  2a2 


where 
R 


JL-    J  ^  j  (ax)     -i 


>C7 


+  4y  Ho(ax) 

a      /o 


x  ni 

dx    +    — x 

2aT 


} 


(17) 


yC? 


Jq  (a(x+z)) 


£7 


dx 


2a' 


2a^ 


H  (a(x+z)) 


ni 


— dx     -    — *■ 

£7  *? 


Jo(a(z-x)) 


J(a(x-z))     -f—  dx     +     -^ 

vE?        2a 


/ 


HQ(a(z-x)) 


>C? 


£7 


dx 


dx 
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We  now  differentiate  the  following  with  respect  to  a: 


q(a)   -  £  Eo(«) 


n 

2^ 


We  get 


2a 


2a 


q'(a)  .  +^  >C7  -  ^    . 

Adding,  we  get 

2  Vn 


(A  +  B  +  C)  +  q'(a)  =  R  . 


When  z  is  not  in  the  neighborhood  of  0  or  1,  R(zja)  vanishes  at  least 

-5/2 

to  order  a     '       as  a  ->oo.  Therefore 


a 

r 


(18)    T(z;a) 


R(zjp)  dp 


oo 


exists.    We  see  then  that  the  left-hand  side  of  equation  (16a)   vanishes  to 
higher  order  than  a     '    as  a  -*co. 

We  now  give  the  value  of  1^  in  closed  form.     We  have 

1 


dK 
3c 


l 
2 


)l-w   I  iawx   iaw( z+x) ,  iawx   iaw I z-x  I 
1  <e    ••  e       +  e    —  e 


dw  • 


J 


We  have  already  shown  that 


lim   G(x,yja)  =  -  |  log  I  *^L 


a  -*•  oo 


Magnus 


W 


has  shown  that  G(x,yja)  is  also  given  by 
1 

>£*   J  eiaw|x-y|  _  eiaw|x+y|  I  dw  # 


i 

2   | 
•o 


) 
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It  follows  that 


...     8K 
a  -*ao 


\  log 


2   2 
z  -  x 

2— 


Prom  the  above  calculations  we  now  get 


(19)      T(zja) 


i 
2 


} 


log 


2   2 

z  -  x 

T~ 


E1(x)  dx  . 


Now  T(Oja)  ■  0, and  so  we  can  write  the  above  as  follows: 

} 

(20)   T(zja)  -  T(0ja) 


2   2 


1 
r 


log|z'-  xt|E1(x)dx  -  -£ 


log  x  E..(x)  dx 


If  we  find  a  solution  for 


(21) 


T(zja)  -  -|    log  |z-x*|  E1(x)  dx 


2  _2 


which  is  valid  for  0  <  z  <  1  then  this  solution  will  also  satisfy  (20). 

2  2 

We  obtain  a  solution  of  (21)  as  follows.  Let  z  -  s  and  let  x  -  t. 


Equation  (21)  becomes 

(21a)  T(  ys"  ja)     =     -  §  log|s-t|  E1(  /F  ) 

Let 

E-^yn 


dt 

2  yt 


2  /r 


Kt)   , 


and  we  get 
(22) 


2iT(  "/s"  ja)     -  log|s-tl  Y(t)  dt     . 
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This  last  equation  has  been  studied  by  T.  Carleman  (cf.  Schmeidler '-'-I ,  p.  56). 
The  solution  is  given  as  follows 

2  /T        2n  log  2  YC  fiZi         I   /s(l-s) 
(23) 

where  the  last  integral  is  a  principle  value  integral. 
We  find  from  (23): 


(210 


E     (v)         1 

2 

T^z>      dz 
?     jo      £7 

Ei^\x.j 

Aog  2  /l-3 

*•»>£?  j 

^     olT(z)  ^ 
?. 

Z    -  X 

o 

(2>Ax)    - 

dE1(x) 
"  dx               * 

dz  , 


(25) 


We  state  this  last  as 


Theorem  3.  Let  the  function  Gj^x)  be  defined  by  (17),  (18), 
(210,  (25).  Then  &,(x)  is  the  first-order  term  arising  from 
a  perturbation  of  the  integral  equation  (?)  for  &>  (x)  at  a  =  oo. 


5,  A  related  result 

We  give  finally  a  result  which,  although  it  is  not  directly  connected 
with  the  diffraction  problem,  follows  immediately  from  it. 

One  of  the  integrals  which  arose  in  the  above  calculation  was 
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sin  aw(x+z)  dw  dx     . 


We  shall  evaluate  this  integral  as  a  — >oo  (in  a  manner  different  from  that 
which  has  preceded) .     We  have 
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If  we  take  the  limit  of  this  last  expression  as  a  -»oo  we  get 
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cos  tz  dt  + 
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sin  tz  dt1 


On  the  other  hand  we  have 
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Equating  these  two  results  we  get 

/"  V*) 

w-    cos  tz  dt     + 

Jo         *^ 


— ^-      sin  tz  dt     =    r     • 

t2  « 

o 


Since  the  sine  transform  of  J,  (t)/t     is  known  we  have  a  formula  which  gives 

2 

the  cosine  transform  of  H.(t)/t  . 

To  our  knowledge  this  has  not  been  evaluated  previously.  We  remark 
that  there  is  some  generality  in  the  method  and  that  other  transforms  can  be 
obtained  similarly.  For  example  if  fl-x      is  replaced  by  (l-x  )  '  ,  etc., we 


can  obtain  in  this  manner 

r 


H2(t) 


cos  tz  dt  + 


J2(t) 


sin  tz  dt1 


h   P(3) 
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